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Abstract. - Optical vortex traps created from helical modes of light can drive fluid-borne 
colloidal particles in circular trajectories. Concentric circulating rings of particles formed by 
coaxial optical vortices form a microscopic Couette cell, in which the amount of hydrodynamic 
drag experienced by the spheres depends on the relative sense of the rings' circulation. Tracking 
the particles' motions makes possible measurements of the hydrodynamic coupling between the 
circular particle trains and addresses recently proposed hydrodynamic instabilities for collective 
colloidal motions on optical vortices. 



A beam of light with helical wavefronts [1] focuses to a ring-like optical trap known as 
an optical vortex [11,14,24], which not only traps micrometer-scale objects, but also exerts 
torque on them [10,13,25]. Since their introduction nearly a decade ago, optical vortices have 
been used to probe the nature of photon orbital angular momentum [5, 10, 19, 20, 25] and to 
create microoptomechanical devices such as rotary pumps and mixers [6,16]. In this paper, we 
demonstrate a new class of microoptomechanical devices resembling Couette cell rheometers 
that are based on optimally matched pairs of concentric optical vortices [12] created with the 
holographic optical tweezer technique [6-8, 18,23]. These paired- vortex machines constitute a 
useful model system for studying many-body colloidal hydrodynamics. 

The field, ip(r) — u(r) exp(i£8), in an optical vortex is characterized by a transverse 
phase profile ipe(r) — Id, where r = (r, 9) is a polar coordinate in a plane normal to the 
optical axis and I is an integer winding number, also known as a topological charge, that 
characterizes the beam's helicity. The amplitude profile u(r) typically is that of a conventional 
TEMoo laser beam whose wavefronts are imprinted with <pe(v) to create the helical mode. In 
our implementation, a liquid crystal spatial light modulator (Hamamatsu X7550 PAL-SLM) 
encodes computer-generated holograms on the wavefronts of a A = 532 nm beam from a 
frequency-doubled diode-pumped Nd:YV04 laser (Coherent Verdi). The modulated beam is 
relayed to a microscope objective lens (Zeiss 100 x NA 1.4 oil immersion S-Plan Apochromat) 
mounted in an inverted optical microscope (Zeiss S100-TV), which focuses it into an optical 
vortex. 
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Fig. 1 - Creating concentric vortices from a beam of light, (a) Gray-scale representation of the phase 
hologram, ip(r), encoding two concetric vortices of topologicl charges h — —30 and h = 80. (b) 
Focused image of the vortices projected by <p(r), after aberration correction, (c) Bright-field image of 
2a = 0.99 /urn diameter polystyrene spheres trapped and rotating in optimized optical vortices. The 
scale bar indicates Wfim. (d) Measured brightness of the laser light around the circumference of the 
rings in (b). (e) Circumferential speeds of an individual colloidal sphere, measured separately on the 
inner and outer rings, and normalized by laser power. 



The same objective lens can be used to create images of objects interacting with the 
projected optical vortex. Provided that care is taken to minimize or correct for aberrations 
in the optical train, this system is capable of trapping one or more colloidal particles and 
circulating them at up to several hundred revolutions per minute. 

If the radial amplitude profile u(r) varies slowly across the optical train's aperture, A, then 
the focused optical vortex takes the form of a ring of radius [5, 26] 

surrounded by a hierarchy of concentric diffraction rings. Here, / is the objective lens' focal 
length and Iq w 10 arises from the SLM's pixellated structure [5]. 

Very recently, Guo et al. [12] showed that the extraneous diffraction rings can be ascribed to 
rays of light emanating from the central region of the phase mask, <pi (r) , while the principal 
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ring is projected from the outer region. This explains why optical vortices with winding 
numbers exceeding £ = 200 can be projected with a pixellated SLM [5, 6, 26] even though 
features in the associated phase mask exceed the device's Nyquist wavenumber near the optical 
axis. It also suggests that the generally undesirable diffraction rings can be suppressed by 
imposing an annular radial profile u(r) = 0(r — ri)9(ro — r) on the source beam, where 
0(x) is the Heaviside step function. The outer limit, r\, replaces the system's aperture A in 
Eq. m to set the optical vortex's radius, and the inner limit, ro, can be adjusted to optimize 
the optical vortex's radial profile. Following Ref. [12], the optimal inner limit that eliminates 
diffraction rings without degrading the vortex's principal ring is approximately given by 

rt^Y^rn, (2) 

where j' e n is the n-th zero of the derivative of the Bessel function Ji{x). Guo et al. [12] point 
out that the principal ring's intensity is diminished when ro exceeds rfj, although its radius 
remains unchanged. 

The unused central region of an optimized vortex's phase mask can be filled with a phase 
profile encoding one or more additional optimized optical vortices. In this case, the inner limit 
of the first phase mask establishes the maximum outer limit of the next, and thereby helps 
to determine the projected radius of the second vortex through Eq. Q). Depending on its 
winding number £, the secondary optical vortex can be either larger or smaller than the first. 
The optimal inner limit of the secondary vortex's annular phase mask again is determined by 
Eq. ©. The number of vortices that can be projected in this way is limited by the size and 
resolution of the SLM. 

A typical phase mask encoding two concentric optical vortices is shown in Fig.^a), and the 
associated pattern of traps appears in Fig. ^b). The two optical vortices' relative intensities 
can be tuned by adjusting their annular phase masks' inner limits. Their relative radii can 
be adjusted by setting both the outer limits and the topological charges. The combination of 
relative intensity, geometry and topological charge helps to establish the relative torques that 
the concentric focused rings of light exert on trapped objects. The resulting tunability is very 
helpful for controlled studies of driven colloidal hydrodynamics. Here we have projected two 
concentric vortices of topological charges l\ = ±30 and £2 = 80, with focused principal rings 
formed at R\ = 2.9 /im and R2 — 6.4 /zm respectively. 

In practice, the part of the input beam passing through featureless regions of the composite 
phase mask in Fig.^a) propagates along the optical axis and focuses to a conventional optical 
tweezer in the middle of the field of view. We have eliminated this undiffracted spot with a 
spatial filter [15] and displaced the concentric vortices by adding a phase function 

^ k (r)=k-r+^p (3) 

where the in-plane wavenumbers, k x and k y , set the in-plane displacement and the axial 
wavenumber k z = 2irz/{\f) displaces the focal plane along the optical axis [6,18]. FigurcQJa) 
shows (/?(r) without fk(r) for clarity. When combined with appropriate control and correction 
of aberrations, these displacements yield the comparatively circular and uniformly bright 
rings in Fig. ^b) . The rings were imaged at low laser intensity by placing a mirror in the 
microscope's focal plane and projecting the reflected light along the microscope's imaging 
train to a charge-coupled device (CCD) camera (NEC TI-324AII). Each ring can trap and 
circulate a single colloidal particle, and the pair can organize multiple particles into concentric 
circulating rings, as shown in the bright-field image in Fig. ^c). 
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Our samples consist of monodisperse polystyrene sulfate microspheres (Bangs Laboratories 
L030305C) 2a — 0.99 in diameter dispersed in a layer of water 18 ± 2/xm thick sandwiched 
between the parallel glass surfaces of a microscope slide and a #1 cover slip. The concentric 
optical vortices' axial intensity gradients were not strong enough to overcome radiation pres- 
sure and trap these particles stably in the axial direction. Consequently, we focused the traps 
near the upper surface, which prevented the particles from escaping while permitting them to 
circulate freely around the rings. 

The concentric vortices' circumferential intensity variations were estimated by fitting the 
radial intensity profile in images such as Fig.^b) to the predicted cross-section [26] . Azimuthally- 
resolved intensity variations, plotted in Fig. ^d), are smooth enough that particles do not 
become trapped in localized hot spots [5] . We measured the individual optical vortices' trans- 
port properties by loading a single particle into one of the rings and tracking its motions [4] . 
Consistent results were obtained for laser powers ranging from P — 0.4 W to P = 2.6 W, 
with the measured single-particle speeds scaling nearly linearly with power. The data from all 
runs were normalized by power and combined into the results plotted in Fig.^e). Somewhat 
paradoxically, regions of minimum circulation speed appear to be correlated with the brightest 
regions of the ring, where the optically-induced torque should be greatest. This demonstrates 
that the intensity variations around each ring's circumference result in attractive optical gra- 
dient forces. The importance of this gradient attraction relative to the torque exerted by the 
photon's orbital angular momentum confirms that the particles absorb a very small fraction 
of the incident photons. 

Despite the particles' preference for brighter regions, there is no evidence of circumferential 
trapping, and we are justified in averaging the circulation speed over angles to characterize 
the rings' overall performance. The mean circulation speed for a single sphere on the outer 
and inner rings is plotted as a function of laser power in Figs. Ufa) and[2Ib). Departures 
from linearity in the power dependence are comparable to our measurement error over the 
range, and can be ascribed to increased hydrodynamic coupling to the upper glass surface with 
increasing radiation pressure. Consistent results for both rings were obtained for l\ = +30 
and i x = -30. 

The hydrodynamic drag on a single sphere of radius a moving at speed v through an 
unbounded fluid of viscosity rj is F v = —70 v, where 70 = 6Trr)a is the Stokes drag coefficient. 
Adding particles to one ring reduces the drag on each sphere through the familiar drafting 
effect. For N spheres equally spaced around a ring of radius R 3> a, the modified single-sphere 
drag coefficient is [22] 



-1 -1 




where fty = {2it/N){j — i) is the angular separation between spheres i and j, and r.y = 
R yj2 — 2 cos Oij is their spatial separation. We can use this result as a starting point to 
predict the filled rings' circulation rates based on our single-particle measurements. 

Coupling between the two rings should increase the drag on spheres in counter-rotating 
rings, l\ = —30, and reduce it for co-rotating rings, l\ = +30. This is indeed the case, 
with the effect being considerably more pronounced for the inner ring (Fig. Ufa)) than the 
outer (Fig. Gib)). The predictions of Eq. Q are plotted as dashed lines in these figures. 
Although Eq. |@J| adequately describes the motions of the outer ring of spheres, it substantially 
underestimates the hydrodynamic coupling among particles on the inner ring. This might 
result from the nonuniform separation between spheres as they circumnavigate the inner 
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Fig. 2 - Colloidal particles driven by concentric optical vortices. (a,b) Rotation rates of the inner 
(Ri) and outer (-R2) optical vortex: single particle in the system (circles), filled vortices co-rotating 
(squares) and counter-rotating (triangles), (c) Hydrodynamic coupling of two concentric, rotating 
rings of particles, measured by the change in the inner ring rotation rate as it's direction is switched, 
Awi. Drag force represented by the rotation rate of the outter ring u)i. Dashed line is the linear fit 
from Eq. JHJ. 



ring's intensity variations, which are more pronounced than on the outer ring, Fig. f^d) and 
(e). 

Even this level of agreement may be coincidental, however. The spheres in this experiment 
are separated from a nearby wall by a center-to-surface distance h that is not accounted for 
by Eq. 0. The no-slip boundary condition on the wall modifies the spheres' far-field flow 
patterns, reducing their mutual hydrodynamic coupling [9]. At the level of the stokeslet 
approximation [21], the equivalent wall-corrected result to O ((/1/i?) 3 ) is 



lN{ R'h 



70 Kh 



8cos( 



N 



l , 3 a Ij (f ) 

. 8R ^ u 



(2-2cos6»i J ) 3 / 2 R 2 



(1 + cos#ij)(l — 3cos#ij) h 
^2 - 2COS01.,-) R 
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is the single particle's wall-corrected Stokes drag coefficient. This still shows an overall re- 
duction in the drag due to drafting, but to a substantially reduced degree. The modified 
functional dependence on the inter-particle separation suggests that the spectrum of insta- 
bilities for such symmetric configurations [22] also will be modified by coupling to bounding 
surfaces. A fit of Eq. (J5J to the results for the outer ring, plotted as a solid curve in Fig.^b), 
yields h = 3a. This is consistent with the measured sphere-wall separation. 

Even at this low level of approximation, the predicted reduction in drag accounts for the 
increased circulation rates when five particles are loaded onto the inner ring and twelve onto 
the outer, as shown in Fig.^c). The circulation rates on the outer ring are increased by 30% 
(Fig. Gib)), and on the inner ring by 47% when the two rings are co- rotating and by 22% 
when they are counter-rotating (Fig. Ufa)), over the range of laser powers applied. 

The data in Fig. [3 also reveal that the many-body hydrodynamic coupling in this system 
depends on whether the two rings of spheres are co-rotating {l\ = +30) or counter-rotating 
(£i = —30). The inner ring, in particular, circulates 20% faster in the co-rotating configura- 
tion. The difference Au>\ in the inner ring's circulation rate increases with laser power, but is 
unlikely to result from optical interactions directly. Rather, the trend can be ascribed to the 
increasing circulation rate of the outer ring with increasing laser power, as shown in Fig.|2c). 
Parameterizing Au>i by lo 2 is reasonable because the outer ring's circulation rate depends only 
weakly on the relative circulation direction. 

We may capture this behavior semi-quantitatively by ignoring the concentric rings' detailed 
structure and treating them instead as concentric cylinders in a Couette geometry. In this 
case, the drag-induced torque between the cylinders is [17] 

T = ^vH-&^(u 2 -^l (7) 

where H is the effective height of the cylinders. The change in frequency upon switching the 
relative circulation direction is then 



3a/ R{ 
2H V'M 



(8) 



Fitting to the data in Fig. Elyields H = 0.15 a. For continuous cylinders with no-slip boundary 
conditions, we would expect H = 2a. The difference can be ascribed to the rings' roughly 25 
percent filling factors, which establish partially sliding boundary conditions. 

This simple model's success suggests that the detailed distribution of particles on the rings 
plays a minor role in establishing the inter-ring hydrodynamic drag, although the discrete- 
ness itself is very important. Consequently, colloidal particles in concentric optical vortices 
should provide a model experimental system for studying how surface roughness influences 
hydrodynamic boundary slip [2,3], with the degree of surface structure being determined by 
the particle size and filling factors. 

We have demonstrated that optimized optical vortices can be used to trap and circulate 
fluid-borne colloidal particles. Concentric circulating rings of particles constitute a model 
test-bed for studying many-body hydrodynamic coupling in mesoscopic systems. Improve- 
ments in the optimized optical vortices' uniformity will make possible detailed investigations 
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of hydrodynamic instabilities in driven many-particle systems. More immediately, the system 
of optically driven colloidal rings shows promise as a microoptomechanical Couette shear cell, 
with potential applications in lab-on-a-chip systems. 

* * * 

This work was supported by the National Science Foundation through Grant Number 
DMR-0451589 and Grant Number DBI-0233971. 
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